The /-electron physics of rare-earth iron pnictides: influence of d-electron 
antiferromagnetic order on the heavy fermion phase diagram 
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Some of the high iron pnictides contain magnetic rare-earth elements, raising the question 
of how the existence and tunability of a d-electron antiferromagnetic order influences the heavy 
fermion behavior of the /-moments. With CeOFeP and CeOFeAs in mind as prototypes, we derive 
an extended Anderson lattice model appropriate for these quaternary systems. We show that the 
Kondo screening of the /-moments are efficiently suppressed by the d-electron ordering. We also 
argue that, inside the d-electron ordered state (as in CeOFeAs), the /-moments provide a rare 
realization of a quantum frustrated magnet with competing J1-J2-J3 interactions in an effective 
square lattice. Implications for the heavy fermion physics in broader contexts are also discussed. 
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PACS numbers: 71.10.Hf, 71.27.-|-a, 74.70.Tx, 75.10.-b 

The homologous rare-earth iron arsenides exhibit anti- 
ferromagnetic (AF) ground states in addition to the high 
temperature superconductivity The sys- 

tems of interest here are the arsenides i?Oa;Fi_a;FeAs, 
with R = Ce, Sm, Nd, Pr,... being magnetic rare 
earths, which have superconducting transition temper- 
atures higher 0,M 0, 0] than the maximal Tc w 26 K of 
La02,Fi_2;FeAs The parent compounds of these sys- 
tems, ROFeAs, have a layered structure, with FeAs and 
RO layers sandwiching each other. They typically show 
a collinear AF order and a structure distortion, which 
are successively suppressed by carrier doping in favor of 
superconductivity [6|. Also of interest are the iron phos- 
phides. LaOFeP was the iron pnictide reported to show 
superconductivity below Tc « 4 K Q. This compound 
has the same layered structure as LaOFcAs, but docs not 
order magnetically [1]. 

The distinction between the iron phosphides and ar- 
senides becomes even more pronounced when La is re- 
placed by Ce. CeOFeP is neither superconducting nor 
magnetically ordered, and its Ce /-electrons exhibit 
heavy fermion behavior with a Kondo temperature Tk ~ 
10 K 01. CeOFeAs has the d-electron collinear AF or- 
130 K. Its /-electrons display a no- 



dering below 



ticeable AF order below T^' « 4 K 0, [T3| , but does not 
show any heavy fermion features. What underlies the 
heavy fermion behavior in CeOFeP and its absence in 
CeOFcAs? One possibility is that this primarily reflects 
the very different interlaycr 3c?-4/ couplings between Ce- 
OFeP and CeOFeAs, as suggested by a first-principle 
LDA-I-DMFT study However, a more complete the- 
oretical estimate using a full density of states, which 
is strongly peaked away from the Fermi energy, sug- 
gests that the effective Kondo couplings in CeOFeP and 
CeOFeAs may in fact be comparable 1^. Muon-spin- 
relaxation and neutron scattering experiments [ij, [14 1 
may also be interpretted in terms of a sizable Kondo 



coupling in CeOFeAs. 

In this Communication, we discuss the possibility that 
the distinction in the d-electron magnetism between Ce- 
OFeP and CeOFeAs plays an important role in influenc- 
ing their heavy fermion behavior. This mechanism is ex- 
pected to play an especially important role when we con- 
sider not only the end materials CeOFeP and CeOFeAs, 
but also the series CeOFeAsi_a;Pa;, which has been pro- 
posed to realize a continuously varying d-electron^AF or- 
der and the associated quantum critical point 



15| 



Studying the effect of the d-electron AF order on the 
heavy fermion phase diagram not only sheds new light on 
the properties of the iron pnictides, but also represents a 
new twist to the heavy fermion physics in general. Typ- 
ically, AF order in heavy fermion metals is induced by 
the RKKY interactions among the /-moments, and the 
heavy fermion phase diagram involves the competition 
between RKKY and Kondo coupling [l^, 17 1. A tunable 



d-electron AF order adds a new dimension to the heavy 
fermion phase diagram. 

In the following, we will consider this effect within an 
extended Anderson lattice model (ALM) appropriate for 
the stoichiometric i?-llll compounds ROFeX (X=As 
or P). The model incorporates the inter-layer hybridiza- 
tion between pnictogen X p-orbitals and rare earth R 
/-orbitals. We note in passing that the derived model 
takes into account the microscopic crystal structure and 
symmetry of the i?-llll compounds. Given that there 
are many materials of the same ZrCuSiAs-type structure 
[isl , with many of them containing magnetic rare-earth 
elements, we expect that our model will also be germane 
to many such related compounds [l^. 

General considerations. The lattice structure of the in- 
fill compound series is schematically shown in Fig. [T] 
Let Fe-atoms be in the {x, y)-plane with the coordinate 
(r, 0), where, r = {ix, iy), ix and iy are both integers (the 
nearest Fe-Fe distance is set to unity). The coordinates of 
X- and i?-atoms are (fp,r]Zp) and {ff^rjZf) respectively. 



2 



(5» n 



-a 




FIG. 1: (Color online) The lattice structure of R-1111 series. 
The small black and red (connected with solid line) circles 
represent Fe and O ions, respectively, and the big blue/brown 
and dashed grey circles are the R and X ions, respectively. 
The small solid red, solid brown, and doted green squares 
describe the Fe, R, and X plaquettes, respectively. Left panel; 
afe-plane; right panel: ac-plane. The dashed and doted lines 
denote Vpf and Vpd, respectively. 



where fp = (i^ + 1/2, iy + 1/2), f/ = {i^ - 1/2, ij, + 
1/2), rj = e'^f'^+'f^ Zp and Zf are the distances of X- 
and i?- atoms to the Fe-plane. We denote the d-, p-, and 
/-electrons by dl^^f), pi^\rp,riZp), and fa"^\rf,rizf), 
with orbital indices a = d^y, d^z, dy^, d^2_j,2, c?3^2_^2, 
H = Px, Py, Pz, and m = 1, • • • , /. 

The model Hamiltonian. The hybridization part of the 
Hamiltonian is given by Hhyj^rid = Hp^ + Hpf , where 

r 

Hpf = E virhP^^^Hrp, ilZp)F!r\rf,riZf) + h.c] . (2) 



Here we introduce Di°'\r) = J2n = da^\ix,iy) + 

d'f \ix + 1, iy) +d''f\ix,iy + l) + di°'\ix + 1, ij^ + 1) and 



p(m) 



- Am), 



H — Hq + Hd + Hf + Hdf + Hint,d + Hint J 
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i^/ = E.^;"""'^[^^i"^^--'/,'7^/)^^r^(r/,7yz/) + h.c.], 
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pd i-p/ ''"V^p'^^- In momentum K-space 
(in lihe reduced Brillouin zone corresponding to two 
Fe-atoms in the conventional cell with lattice constant 



= V2), 



Hd - Z^K ^d 9d (^)%K<TVKa' 



Z^K "^d/ I'd/ V-^-'l"r,Ko-Jr)'Ko- 
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and 



+ h.c], where d. 



(a) 

and /^Kcr ^'^'^ ^1"^^ Fourier transform of d- and /- electron 
operators in the sublattices rj = A or B, respectively. 
The K-dependence of the dispersions and d-f hy- 
bridization is only encoded in the form factors, given 
by g^/^\K) = g^f'^K) = 4 + 2cos{K,a) cos{Kya), 
g(/^\K) = gf^\K) = 8cos(if,a/2)cos(if,a/2), 

ffrf^ '(^J = Scos-iA^^a/Zjcos^Aya/Zj, ffif^^CK) = 

g^^^^(K) = 8cos(i^^a/2)cos2(A'ya/2). 

The d-electron correlations. For moderate large Ud, 
we may start from the strong coupling limit yielding the 
frustrated J1-J2 Hcisenbcrg model for the d-electrons [l^, 
22| . The itinerancy of the c?-electrons will further 



g/(K) = 16cos2(i^^a/2)cos2(/i:j^a/2), 
SBB)^^^ = 8cos2(ii:^a/2)cos(ii'ya/2), 
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1/2, 772/)-f /rn*x+3/2, z.,+1/2, 77Z/)+/rH*x + l/2, z,- 
1/2, rjZf) + fa"^\ix + 1/2, iy + 3/2, r]Zf) as the plaquette 
operators of d- and /-electrons around X-atoms. (Sum- 
mations over the repeated spin and channel indices are 
implied hereafter unless otherwise specified.) 

The interaction part of the Hamiltonian, Bint = 
Hint,d + Hint,p + Hintj; coutaius the usual on-site 
Coulomb interactions {Up, Ud, and Uj) and the Hund's 
coupling {JH,d)- The total Hamiltonian is then B = 
Bq + Biiyi,rid + Bint 1 with Bq containing the primitive 
site energies of d-, p-, and /-electrons denoted by e^"'. Balm 
Ep^^ , and e^p^ , respectively. 

It is expected that Up is small compared to the other 
Coulomb interactions. We will therefore set Up = 0, in 
which case the p-orbitals can be readily integrated out. 
The obtained effective Hamiltonian B takes the form 



reduce the ordered moments and eventually lead to a 
paramagnetic phase [l5| . In fact, both the weak- and 
strong-coupling limits suggest that the staggered magne- 
tization Md = E„ M^") = -(1/iV) EK({<"^j<1^+QT - 

'^r;Ki'^r;K+Qi}) ^ dominating order parameter with 
Q = (tt, tt) and N being the number of K points in the 
reduced Brillouin zone. For the purpose of demonstrat- 
ing the effect of d-electron order on the Kondo effect, we 
treat M^"'' as the mean field parameters and approxi- 
mate B,nt,d by ^dEKf^^^d"^KKl'^!,K+Q<T + li-c-], with 
Jd being the effective coupling strength. The AF order- 



JdMjf^ , is sizable for FeAs but vanishes 



for FeP. 

Kondo effect vs. d-electron ordering. In order to un- 
derstand the competition between the Kondo effect and 
d-electron AF order, we first neglect the /-electron or- 
dering. We are then led to consider 
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In the absence of d-electron ordering, Eq.(4) is th 
ALM with weak /-electron dispersion and momentum 
dependent hybridization. (The effect of momentum 
dependent hybridization on the Kondo effect has recent! 
been studied in other contexts [1^, HSl-) For sufScientl 
large Uf, and with the /-levels being well below Ferm 
energy, we are in the Kondo limit. 

To concretely demonstrate how the d-electron AF oi 
der influences the Kondo effect, we consider the resultin, 
Kondo lattice model with a single /-electron channel am 
two c?-electron bands. In the slave-boson representatior 
this becomes 

k ^ ^ k ^ 

k 

-ljKY.v,fi\^)bMlA1 +ii.c.]. (5) 

^ k 

Here, the Lagrange multiplier A enforces the single oc- 
cupancy of /-electrons. The mean- field parameter ba = 
{flad^J)/2 describes the Kondo screening and sets the 
Kondo scale, Tk oc b^. The anisotropic hybridization 
form factor Vd/(k) = 4cos/ca:/2cosfcy/2. The energy dis- 
persion for d-electrons are taken to 136 pHj : 





(k) = 


—2ti cos kx — 2t2 cos ky — At^ 


cos kx 


cos ky , 


,(2) 


(k) = 


—2t2 cos kx — 2ti cos ky — At^ 


cos kx 


COS ky , 


g(12) 


(k) = 


e^^"'^-'(k) = — 4t4 sin fca: sin fcy , 






with 


ti = 


^Ij ^2 = 1.3, = = - 


-0.85. 


In our 



numerical study, we choose Jk — 0.04, temperature 
T = 10-i°|ti|, and the lattice size Nl = 3200 x 3200. 
When we vary the AF order parameter, the chemical po- 
tential is adjusted such that the d-electrons are fixed at 
the half-filling Ud = 2.0. 

Fig. [2] shows that the d-electron AF order rapidly sup- 
presses the Kondo scale. This suppression is closely re- 
lated to the depression of the d-electron density of states 
(DOS) in the coUinear AF state of undoped iron arsenides 
(see the inset of Fig. ^ . The feature of low energy DOS is 
sensitive to the degree of nesting and the DOS minimum 
is not necessarily located precisely at the Fermi energy 
(see, e.g., the case of Md = 0.021 in Fig. [2]); the latter 
explains the effective Kondo scale first rising and then 
dropping with the AF order. Furthermore, the incom- 
plete nesting of the Fermi surface keeps the depressed 
DOS finite (unlike, say, in the superconducting state) at 
the Fermi energy such that the T = ground state has 
the /-moment always Kondo screened on the lattice. 

We should stress that, for the purpose of a semi- 
quantitative assessment of the proposed mechanism, we 
have considered the upper limit for the Kondo scale in 
the AF state: wc have coupled the /-moments to only the 




FIG. 2: (Color online) Mean-field Kondo parameter & as a 
function of the Q = (tt, 0) (in the notation of the one-Fe Bril- 
louin zone) staggered magnetization Md. The Kondo temper- 
ature Tk oc b^. Md is measured in /is/Fe. The inset shows 
the d-electron density of states for Md ~ (dotted- black), 
0.021 (black-red), 0.165 (dashed-green) , 0.296 (dash-dotted- 
blue). 



quasiparticles of the d-electron AF state and have also 
neglected the /-moment ordering; moreover, a genuine 
/-electron quantum phase transition will be induced by 
breaking the Kondo screening upon the inclusion of the 
standard RKKY-Kondo competition |26„^,^. We can 
therefore infer that the mechanism proposed here pro- 
vides a viable basis to understand the distinct /-electron 
heavy fermion behaviors in CeOFcP {Md « 0) and Ce- 
OFeAs {Md w 0.8[l3|). Our results also set the stage for 
understanding the evolution of the heavy fermion behav- 
ior in the CeOFeAsi_a;Pa; series. In general, there will 
be two magnetic quantum critical points Xci and , as- 
sociated with the d- and /-electrons, respectively. The 
RKKY interaction would then dominate in the interme- 
diate region of x, leading likely to a ferromagnetic order 
before the heavy fermion state is approached. 

Magnetic frustration of the f -electrons. We now turn 
to the exchange interactions among the /-moments. Con- 
sider first the supcrcxchange interaction, which can be 
derived by integrating out the virtual valence fluctua- 
tions of the /-electrons. From Eq. we end up with 

Hf - Er~jt-'"''^S^^\rf,vzf) ■ s'^'-^rf^rjzf), where 

Sp^\r,r]Zf) = X^n 'S/™'' (rf ' ^^/) ^^'^ summations of /- 
electron spins in the corresponding plaquettes associated 
with f, and jj."'"') « 2[y;"'"')]2( L_ _ This 

is the supcrcxchange interaction associated with the R- 
X-R path, which does not mix the odd and even sublat- 
tices of the /-sites in a single RO layer (see Fig. [Hla)). 
There will also be a superexchange interaction from the 
R-O-R path, due to the hybridization between the 4/- 
orbitals of X-atoms and the 2p-orbitals of 0-atoms; this 
superexchange mixes the odd and even sublattices (see 
Fig. EKb)). In the notations of an effective square lat- 
tice of the /-sites {c.f Fig. [3^c)) the R-O-R path gives 
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(b) 
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FIG. 3: (Color online) The would-be ordering patterns of 
the /-electrons due to the superexchange interactions via (a) 
the R-X-R process alone or (b) the R-O-R process alone, 
(c) illustrates the combined exchanges, viewed as J1-J2-J3 
interactions of an effective square lattice within an RO layer, 
which are expected to turn the orders of (a) and (b) into a 
helical one. The blue and brown circles label the rare-earth 
sites in the same way as in Fig. [T] 



rise to the nearest-neighbor (n.n.) interaction J^'^ and 
the next-nearest-neighbor (n.n.n.) while the R- 

X-R path yields the n.n.n. J2 and the third-nearest- 
neighbor (n.n. n.n.) J^^'' ■ (Note that J^"^'' and J^^^ 
correspond to the n.n. and n.n.n. interactions in 
the odd/even sublattices separately.) The resulting /- 
electron spin Hamiltonian becomes a J1-J2-J3 Heisenberg 
model (Fig.ElJc)) 



(6) 



where Ji = J2 



jf)+jf\and J3 = Jf\ 



In this way, the /-moments of CeOFeAs provides a re- 
alization of a geometrically frustrated quantum magnetic 
system in two dimensions. Quantum frustrated magnets 
have been the subject of theoretical studies for a long 
time, and continue to attract extensive interest [2^. 
However, suitable materials with spin-1/2 are rare. In 
this context, it will be very important to clarify the mag- 
netic behavior of the /-moments in CeOFeAs and related 
arsenides. 

The Sf ■ Sf form given earlier corresponds to 



being equal to 1/2, and further bond-angle 

t{X) I j{X) 



considerations imply that J^^"^' /J2^' will be larger than 
1/2 but still not far away from it. Similar considerations 
would suggest that j'^^ / j'^^ w 1/2. We will therefore 
expect J2 > Ji/2 and a sizable J3/J1. In this range, 
the Neel and coUinear orderings are excluded. Instead, 
an incommensurate helical phase with the ordering vec- 
tor (g, tt) or (g, a) is the most likely ground state, where 
cosg = ^"^Ijf^ [3Q(]. Neutron scattering and muon spin 
relaxation experiments in polycrystal CeOFeAs appear 
to have seen a hehcal /-electron ordering l^, l^ - 

In the d-electron paramagnetic regime, there will also 
be an RKKY interaction. The latter is expected to be 
ferromagnetic given the relatively small size of the Fermi 
surfaces, and this is consistent with the enhanced ferro- 



magnetic fluctuations of the heavy fermion state observed 
in CeOFeP [9|. Still, the frustrating J1-J2-J3 superex- 
change interactions will continue to operate, helping to 
suppress the tendency for AF ordering. 

Discussion and summary. A number of other conse- 
quences of the p-f hybridization are relevant to the iron- 
pnictides phase diagram. First, in the heavy fermion 
phase, the momentum-dependence of the induced d- 
f hybridization will generally smear the hybridization 
gap (which has nodal lines along K^, = ztn/a and 
Kj, = ±7r/a), and this could be visible in the optical- 
conductivity spectrum. Second, the induced d-f hy- 
bridization depends on the Fe-AT and X-R distances. In- 
creasing pressure along the c-axis will decrease the dis- 
tances and increase the hybridizations, and eventually 
enhance Tr- [ii|. Finally, in light of the fact that the /- 
electron ordering is further suppressed by the competing 
J1-J2-J3 interactions, the transition or crossover from the 
superconducting to the heavy fermion phases may take 
place at sufficiently high pressures in the carrier-doped 
superconducting materials. 

In summary, we have considered a mechanism for 
weakening the Kondo screening effect through the an- 
tiferromagnetic order of the conduction electrons, and 
implemented it in an extended Anderson lattice Hamil- 
tonian. For the iron pnictides, our mechanism is 
semi-quantitatively viable to explain the observed ex- 
istence/absence of heavy fermion behavior in CeOFeP 
and CeOFeAs, respectively. More broadly, our mecha- 
nism goes beyond the standard picture of heavy fermion 
physics, viz. the RKKY and Kondo competition, and can 
therefore shed new light on the phase diagram of heavy 
fermion systems in general. Finally, we have proposed 
that the /-electrons in the parent iron arsenides represent 
a rare model system for quantum frustrated magnetism 
in two dimensions. 
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